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Abstract 

In  the  binomial  tree  model,  we  provide  efficient  algorithms  for  computing  an  accurate  lower  bound 
for  the  value  of  a  European-style  Asian  option  with  either  a  fixed  or  a  floating  strike.  These  algo¬ 
rithms  are  inspired  by  the  continuous-timeanalysis  of  Rogers  and  Shi  [10].  Specifically  we  consider 
lower  bounds  on  the  option  value  that  are  given  by  the  expectation  of  the  conditional  expectation  of 
the  payoff  conditioned  on  some  random  variable  Z.  For  a  specific  Z,  Rogers  and  Shi  estimate  this 
conditional  expectation  numerically  in  continuous  time,  and  show  experimentally  that  their  lower 
bound  is  very  accurate.  We  consider  a  modified  random  variable  Z  that  gives  a  strictly  better  lower 
bound.  In  addition,  we  show  that  this  lower  bound  can  be  computed  exactly  in  the  n-step  binomial 
tree  model  in  time  proportional  to  ri’ .  We  show  that  computing  the  approximation  is  equivalent  to 
counting  paths  of  various  types,  and  that  this  can  be  done  efficiently  by  a  dynamic  programming 
technique.  We  present  other  choices  of  Z  that  yield  accurate  and  efficiently-computable  lower 
bounds.  We  also  show  algorithms  to  compute  a  bound  on  the  error  of  these  approximations,  so  that 
we  can  compute  an  upper  bound  on  the  option  value  as  well. 
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1  Introduction 


The  binomial  tree  model  of  Cox,  Ross  and  Rubinstein  [3]  is  commonly  used  as  a  discrete-time 
approximation  to  a  continuous-time  geometric  Brownian  motion.  The  pricing  of  simple  calls  and 
puts  in  this  model  can  be  done  efficiently.  However,  path-dependent  options  such  as  Asian  options 
are  notoriously  hard  to  price  in  this  model.  The  Black-Scholes  differential  equation  for  the  price  of 
such  options  has  no  known  closed  form  solution. 

The  payoff  of  an  Asian  option  is  based  on  the  arithmetic  average  of  the  stock  price.  In  particular, 
suppose  the  option  expires  at  time  T  and  let  Aj  denote  the  arithmetic  average  of  the  stock  prices 
up  to  time  T.  Then  the  payoff  of  a  European  style  Asian  call  with  strike  L  is  (.47  —  L)'^  where  a:"*" 
denotes  max{a;,  0}.  Similarly  thepayoffof  aEuropean-styleAsian  put  with  strike  L  is  (L  —  .47  )+. 
For  the  so-called  “floating-strike”  versions  of  these  options,  the  L  is  replaced  by  Sr,  the  stock  price 
at  time  T.  Such  options  are  of  obvious  appeal  to  a  company  which  must  buy  a  commodity  at  a  fixed 
time  each  year,  yet  has  to  sell  it  regularly  throughout  the  year  [13].  These  options  allow  investors  to 
eliminate  losses  from  movements  in  an  underlying  asset  without  the  need  for  continuous  rehedging. 
Asian  options  are  commonly  used  for  currencies  [13],  interest-rates  and  commodities  such  as  crude 
oil  [5]. 

Several  researchers  have  devised  approximation  algorithms  for  Asian  options.  This  work  has 
been  of  three  general  types.  Monte  Carlo  methods  and  Fourier  transform  techniques  have  been 
explored  by  Kemna  and  Vorst  [6],  and  Carverhill  and  Clewlow  [2]  respectively.  The  second  category 
of  approaches  consists  of  attempts  to  approximate  the  distribution  of  the  arithmetic  average  (for 
which  no  simple  specification  is  known)  by  a  more  tractable  one.  Examples  of  such  work  include 
those  by  Ruttiens  [11],  Levy  [7],  [8],  Levy  and  Tlimbull  [9],  and  Thmbull  and  Wakeman  [12].  The 
third  approach,  considered  by  Yor  [14]  and  Geman  and  Yor  [4],  is  to  derive  formulas  for  theLaplace 
transform  of  an  Asian  option.  However,  there  have  been  no  numerical  studies  of  the  inversion  of  this 
Laplace  transform,  and  no  simple  analytical  inversion  has  been  found.  None  of  the  above  methods 
offers  an  approximation  for  the  Asian  option  price  that  is  both  fast  and  very  accurate. 

The  work  in  the  present  paper  is  most  closely  related  to  that  of  Rogers  and  Shi  [10].  They 
provide  a  method  for  computing  lower  bounds  on  the  price  of  an  Asian  option  in  the  continuous¬ 
time  geometric  Brownian  motion  model.  Their  main  idea  is  the  following.  Let  be  the  final 
payoff  of  a  European  style  Asian  (call  or  put)  option.  For  instance,  for  a  European  style  Asian  call 
option,  Wr  =  .47  —  L.  Then  for  any  random  variable  Z,  we  can  lower-bound  the  option  value 
using  Jensen’s  inequality: 

E(iy+)=E(E(iy+lZ)) 

>E(E(iy7|^)+)  (1) 

The  lower-bound  on  the  right  in  (1)  can  be  used  as  an  approximation  for  E(Wj^).  The  quality  of 
the  approximation  clearly  depends  on  the  choice  of  the  random  variable  Z.  This  approximation  is 
an  integral  in  continuous  time,  and  Rogers  and  Shi  estimate  this  integral  numerically  for  various 
choices  of  Z.  In  the  present  paper,  we  show  that  for  several  choices  of  the  random  variable  Z  we 
can  efficiently  compute  the  approximation  (1)  exactly  in  the  binomial  tree  model.  We  show  this  for 
random  variables  Z  that  are  related  to,  but  not  exactly  the  same  as,  the  ones  considered  by  Rogers 
and  Shi. 
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One  way  to  think  of  the  above  approximation  in  the  w-step  Binomial  tree  model  is  the  following. 
We  let  W,i  and  denote  the  random-variables  corresponding  to  the  continuous-time  random- 
variable  W-r  and  A'l  respectively.  The  difficulty  in  computing  the  exact  option  value  E{W,f ) 
comes  from  the  fact  that  nearly  every  stock  price  path  has  a  different  (arithmetic)  average  stock 
price  (and  therefore  a  different  and  there  are  2"  paths.  The  approximation  (1)  can  be  thought 
of  as  the  value  of  an  option  with  Wn  replaced  by  W,[,  defined  as  follows.  Let  C,-  be  the  collection 
of  all  paths  with  the  same  value  of  the  random  variable  Z  =  Zi.  Now  let  the  W',^  on  every  path 
in  Ci  be  equal  to  the  average  value  of  over  the  collection  C,.  In  other  words,  in  using  (1)  to 
estimate  the  option  value,  we  are  making  the  approximation  that  all  paths  with  Z  —  Zi  have  the 
same  average  stock  price,  equal  to  the  average  of  the  average  stock  price  over  these  paths.  This 
view  will  he  useful  to  understand  intuitively  why  some  choices  of  Z  yield  better  lower  bounds  than 
others. 

Rogers  and  Shi  find  that  a  particularly  good  lower  bound  is  obtained  by  taking 

Z  =■  (  Bu  du. 

Jo 

which  is  T  times  the  average  value  of  the  Brownian  from  time  0  to  time  T.  In  the  n-step  binomial 
tree  model,  the  Brownian  average  corresponds  to  the  average  of  the  underlying  random  walk.  The 
steps  of  this  random  walk  are  given  by  the  random  variables  X,  =  ±1,  i  =  1, 2, . . . ,  ra.  Ihe 
position  of  this  random  walk  at  time  k  is  given  by  Iq  =  0.  and 

n  =  fc  =  l,2,...,n. 

i=l 

The  average  position  from  time  0  to  time  n  is  Znjn  where 

k 

Zk  =  YlYi,  k  =  0,2,...,n. 
i=0 

The  problem  with  choosing  Z  =  Z,fn  is  that  one  cannot  compute  the  expression  (1)  efficiently 
in  the  binomial  model.  However,  if  we  take  Z  =  (Z„,  S„)  where  5,,  is  the  stock  price  at  time  n, 
then  (1)  is  a  strictly  better  approximation,  and  it  can  be  computed  exactly  in  time  proportional  to  nJ 
(Section  5).  We  also  show  by  means  of  numerical  experiments  (Section  7)  that  the  approximation 
(1)  for  this  choice  of  Z  agrees  excellently  with  the  exact  option  value  in  the  binomial  model. 

Note  that  the  geometric  stock  price  average  is  a  simple  function  of  Z„ .  This  helps  explain  intu¬ 
itively  why  conditioning  on  {Z„,  S„)  yields  a  good  approximation:  paths  with  the  same  geometric 
stock-price  average  that  end  at  the  same  stock  price  will  likely  have  arithmetic  stock-price  averages 
that  do  not  differ  too  much.  The  reason  why  this  approximation  can  be  computed  in  time  propor¬ 
tional  to  only  (rather  than  2'*)  is  that  the  number  of  possible  values  of  is  only  proportional 
to  n^.  By  contrast  the  number  of  possible  values  of  the  stock  price  average  .4,,  is  2".  As  we  show 
in  Section  5,  this  allows  use  to  use  dynamic  programming  to  compute  the  approximation  in  time 
proportional  to  nJ. 

In  this  paper  we  also  consider  Z  =  S^.  8^.  •  •  •  ,  ■S'nd<  Sn,  where  0  <  ni  <  n-2  <  ■  ■  ■  <  Ud  < 
n  are  fixed  integers.  We  are  thus  conditioning  on  the  stock  prices  at  certain  fixed  times.  We  show 
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(Section  4)  that  for  a  fixed  d,  the  expression  (1)  can  be  computed  exactly  in  the  binomial  tree  model 
in  time  proportional  to  Thus  as  we  increase  d  we  get  increasingly  better  approximations,  if 

we  are  wilHng  to  spend  more  computation  time.  In  Section  6  we  describe  a  hybrid  approximation 
where  we  condition  on  Zn  in  addition  to  these  stock  prices.  We  show  that  for  all  our  choices  of  Z, 
our  algorithms  require  time  proportional  to  some  constant  power  of  n.  Thus  our  running  times  are 
polynomial  in  n  rather  than  exponential. 

Rogers  and  Shi  show  that  the  error  in  the  above  approximation  is  bounded  above  by 

iE[var{IFlZ)^].  (2) 

We  show  that  this  bound  can  also  be  computed  efficiently  for  our  choices  of  Z.  Thus  we  can 
efficiently  compute  both  a  lower-bound  and  an  upper  bound  on  the  option  price.  Rogers  and  Shi  do 
not  show  how  to  compute  their  approximations  or  errors  efficiently  in  the  binomial  tree  model. 

The  paper  is  organized  as  follows.  In  Section  2  we  describe  the  binomial  model  and  related 
notation.  In  Sections  3,  4,  5  and  6  we  describe  polynomial-time  algorithms  to  compute  the  lower 
and  upper  bounds  on  the  price  of  an  Asian  option.  Throughout  we  illustrate  our  ideas  for  a  fixed 
strike  Asian  call,  and  it  should  be  clear  how  to  extend  these  to  Asian  puts,  and  floating-strike  Asian 
options.  In  Section  7  we  show  numerical  results. 

2  The  binomial  model  for  stock  prices 

In  this  paper  we  consider  an  Asian  option  written  on  an  underlying  risky  asset,  which  we  refer  to  as 
the  “stock”.  Following  Black  and  Scholes  [1],  we  make  the  standard  assumption  that  the  price  of 
the  stock  at  time  t  satisfies  the  stochastic  differential  equation  for  geometric  Brownian  motion: 

dS(t)  =  dt  +  aS{t)  dB{t), 

where  the  drift  rate  p,  and  volatility  cr  are  constant.  B(t)  is  a  Brownian  motion  process.  We  will 
also  assume  that  the  risk-free  interest  rate  is  a  constant.  The  derivative  security  expires  at  time  T. 

The  continuous-time  function  S  (t)  can  be  approximated  in  the  form  of  a  Cox-Ross-Rubinstein 
[3]  binomial  tree,  as  follows.  The  life  of  the  option  is  divided  into  n  time  steps  of  length  £\t  =  T /n. 
In  time  At  the  stock  price  moves  up  by  a  factor  u  with  probability  p,  or  down  by  a  factor  d  =  lju 
with  probability  q  =  \  —  p,  where 

u  —  e  ^  , 

—  d 

P= - 

u  —  d 

Time  k  in  this  discrete  model  corresponds  to  continuous-time  fcAt.  Also  we  assume  that  a  unit  of 
currency  lent  or  borrowed  at  time  k  will  be  worth  1  -I-  r  units  at  time  k  +  1.  For  brevity  we  write  R 
to  denote  1  -|-  7’.  This  binomial  tree  model  will  be  assumed  throughout  the  paper. 

More  formally,  we  have  a  sample  space  of  all  possible  sequences  of  n  coin-tosses  (i.e.  H’s 
and  T’s),  where  an  H  corresponds  to  an  up-tick  of  the  stock  and  a  T  corresponds  to  a  down-tick.  A 
typical  sequence  (or  “path”)  in  Q  is  written  ui  =  . . . w„,  where  represents  the  ith  coin-toss 
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and  u>i  6  {H,T].  A  /:-prefix  of  a  path  w  is  the  subsequence  u>iu>2  ■  ■  •  w*,.  constisting  of  the  first  k 
coin-tosses.  For  =  1,2, . . . ,  n,  define  the  random  variable  X*  on  fl  as  follows: 


We  will  need  the  following  random  variables: 

Yo  ==  0, 

j=l 

i=0 


if  W/t  =  H, 

if  Ljk  —  'Y' 

k  >  1, 

(3) 

k  >  0, 

(4) 

Also  let  Hk  (w)  denote  the  number  of  heads  in  the  fc-prefix  of  w,  i.e.,  the  number  of  heads  up  to 
and  including  time  k.  We  define  Ho(oj)  =  0  for  all  w.  In  this  paper  it  will  be  useful  to  consider  the 
binomial  lattice  where  a  point  with  coordinates  {k,  h)  represents  all  paths  w  with  =  h.  We 
say  that  the  value  of  a  random  variable  (e.g.,  stock  price  Sk)  at  the  point  [k,  A)  is  ®  if  the  value  of 
the  random  variable  at  time  k  is  a;  on  a  path  w  such  that  =  h,  i.e.,  the  value  of  the  random 

variable  at  time  k  is  x  when  there  have  been  k  heads  up  to  and  including  time  k. 

The  stock  price  at  time  k  is  denoted  Sk,  A  =  0, 1, ....  n,  and  is  given  by 


Thus  the  stock  price  at  point  (A,  A)  is  .  The  average  stock  price  at  time  k  is  defined  as 


.4*  —  {So  -b  5i  -t- . . .  -b  Sk}/{k  4- 1) . 

The  payoff  of  an  Asian  call  with  strike  L  at  time  n  is  (.4,,  —  L)+  and  the  payoff  of  an  Asian  put  is 
(L  —  .4,,)+.  The  payoff  of  a  floating-strike  Asian  call  at  time  n  is  (.4„  —  5’,,)+,  and  for  a  floating- 
strike  Asian  put,  the  payoff  is  (S,,  —  .4,,)+.  In  general  we  write  W,f  to  denote  the  payoff  of  any 
of  the  above  options  at  time  n.  The  value  of  such  an  option  (at  time  0)  is  defined  as 

V  =  E(W+)/i?'‘. 

In  this  paper  we  are  concerned  with  estimating  £(!!''+ ),  and  we  will  ignore  the  discounting  factor 
R". 

Let  T  denote  the  cr-algebra  consisting  of  all  subsets  of  fi,  and  Tk  denote  the  cr-algebra  gener¬ 
ated  by  the  first  k  coin-tosses,  i.e.,  by  A'l , . . .  ,Xk-  We  will  always  assume  the  probability  measure 
P  on  (fl,  X)  given  by: 

P[X/t  =  1]  =  p,  A  =  1, 2, ... ,  n. 

Thus  for  any  random  variable  Z,  Y,{Z\Tk)  denotes  the  conditional  expectation  of  Z  given  the  fist 
k  coin  tosses.  A  random  variable  X  is  .7^/; -measurable  if  it  depends  only  on  the  first  k  coin  tosses. 
An  Xo-measurable  random  variable  is  non-random.  A  sequence  of  random  variables  {Zk}%-o  is 
called  a  process,  and  we  will  often  refer  to  it  simply  as  “the  process  {Zk]".  A  process  {Zk}  is  said 
to  be  adapted  if  each  Zk  is  X/j-measurable.  For  instance  the  stock  price  process  {5/;}  is  adapted. 
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3  Approximations  based  on  conditional  expectations 

In  this  and  the  remaining  sections  we  describe  algorithms  for  approximating  the  value  of  a  European- 
style  Asian  call  option  in  the  binomial  tree  model.  As  mentioned  in  the  introduction,  all  the  algo¬ 
rithms  are  inspired  by  the  continuous-time  analysis  of  Rogers  and  Shi  [10].  In  this  section  we 
describe  a  general  scheme  for  computing  the  approximation  (1)  and  the  error  bound  (2).  In  the  next 
three  sections  we  show  how  to  apply  this  scheme  to  specific  random  variables  Z. 

In  all  of  our  approximations,  the  random  variable  we  condition  on  will  be  a  combination  of 
and  one  or  more  other  random  variables,  which  we  collectively  denote  by  Z.  Our  approximations 
(i.e.,  lower  bounds)  are  therefore  of  the  form 

E{E(Wn\Z,Hn)+),  (5) 

where  the  random  variable  Z  could  be  multi-dimensional.  The  reason  we  always  condition  on 
is  that  all  paths  with  the  same  value  of  have  the  same  probability:  this  aUows  us  to  reduce  the 
computation  of  the  expectation  (5)  to  a  path-counting  problem. 

To  describe  how  to  compute  (5),  we  introduce  some  useful  notation.  Let  N{k,  z)  denote  the 
number  of  paths  with  =  k  and  Z  =  z.  Equivalently,  N (fc,  z)  is  the  number  of  paths  with 
Z  =  z  that  go  through  the  lattice  point  (n,  k) .  Let  Mi  {k.z:  t,  h)  denote  the  number  of  paths  with 
Z  =  z  that  go  through  (n,^•)  and  (Lh).  Let ti.hii  t2, /12)  be  the  number  of  paths  with 
Z  =  z  that  go  through  (n,  k),  (ti.hi)  and  (t2,  /i2)-  Let  C(k}  denote  the  set  of  possible  values  of 
the  random  variable  Z  on  paths  with  H„  =  k.  Then  for  an  Asian  call  our  approximation  can  be 
written 


E(E(I7,.  I 

=  EiE{An-L\Z,H,y) 

=  E(E(.4„ 


=  ^  E(Z=z,H,,  =  k) 

h=QzeC{k) 


=  E  E  pU'‘-'-’N{k,z) 

k=0  3ec{k) 


E(.4„  \Hn  =  k,Z=z)-L 

+ 


E(.4„  I  Hn  =  k,Z  =  z)-L 


=  E  E 


z)E(.4„  \H„  =  k,Z  =  z)-  N(k,  z)L 


(6) 


k=0:iec(k)  '■  ■■ 

If  the  size  of  the  sets  C{k)  is  bounded  by  a  constant  times  /(n),  then  the  number  of  iterations  of 
the  two  summations  above  is  proportional  to  n  f{n).  Note  that 


W{k,  z)  =  N(k,  z)E(.4„  \H„=k,Z  =  z) 


is  the  sum  of  .4,,  over  all  paths  with  H„  =  k  and  Z  =  z.  When  we  write  down  the  expression  for 
this  sum,  we  notice  that  the  stock  price  at  the  lattice  point  (t,  h)  (which  is  appears  in  this 

sum  with  a  coefficient 


Ml  {k.  z;  t,  h) 


(»+l) 
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so  we  can  write  W {k,  z)  as 


'.”  +  <=0  A=o 

Thus  the  approximation  for  the  Asian  call  is  given  by  the  expressions  (6)  and  (7).  To  complete  the 
description  of  the  approximations,  we  only  need  to  show  how  to  compute  the  quantities  N{k,  z), 
Mi(k,  x;  t,  k)  and  the  sets  C{k).  In  the  following  sections  we  do  this  for  three  specific  random 
variables  Z.  Suppose  the  computation  of  each  M^(k,  z\  t,  h)  takes  time  proportional  to  ri{rr). 
Then  since  the  time  to  compute  each  M]  {k,  z;  t.  h)  dominates  that  required  to  compute  N  (fc,  z) ,  it 
is  easy  to  see  that  the  overall  approximation  takes  time  proportional  to 

n^f(n)T,{n).  (8) 

We  now  consider  the  computation  of  the  error  bound  (2).  Again  in  the  case  of  an  Asian  call  we 
can  write  the  expectation  in  the  error  bound  as 

E[var(iy,.  I  Hn,Z)^]  =  E[var(.4„  -  L  |  Hn,Z)^ 

=  E  E  nHn  =  k,Z  =  z)[y^v{An-L\Hn,Z)i] 

*=0  z€C{k) 

*=0  zec{k) 


where 


var{.4„  ~L\H,„Z)  =  E[(.4„  -  Lf  \  Z]  -  [E(.4„  -  L  [  Zf.  (9) 

We  already  know  how  to  compute  the  second  term  in  (9)  above,  and  the  first  term  is 

E[(.4„  -  Lf  I  2]  =  E{Al  ~  2.4„L  +  \  Z] 

=  E{Al  1  H,„  Z)  +  L^-  2LE(4„  |  Z), 

where  we  already  know  how  to  compute  the  last  term,  so  we  only  need  to  show  how  to  compute 
N{k,  z)E{Af^  I  H,^,  Z).  This  is  the  sum  of  .4^,  over  all  paths  with  if,,  =  k  and  Z  =  z.  When  we 
write  down  this  sum,  we  notice  that  there  are  two  kinds  of  terms.  The  first  kind  are  those  involving 
the  square  of  a  stock  price.  The  square  of  the  stock  price  at  lattice  point  (t,  h)  (which  is 
appears  in  this  sum  with  a  coefficient 

t,h) 

(n  +  l)^ 

The  remaining  terms  involve  cross  products  of  stock  prices  at  different  points  on  the  lattice.  The 
product  of  the  stock  price  at  (tj./ii)  and  ((2,  /12),  ti  7^  t2  (which  is  appears 

with  a  coefficient 

2M2(fc,z;  t2.^2) 
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Thus  we  may  write 


*=0A=0 


+  2 


^2  n  ti  n  «2 


(n  +  1)^ 


'll  H  Y  -^^2(fc,z;  h.ki-,  t2,k2)u 


2hi-hy2h2-ti 


^  (10) 


<1=0  /ii=0  (2=«1+1  /12=Ai 

It  is  straightforward  to  check  that  the  time  required  to  compute  the  error  bound  is  proportional 


to 


n^f{n)T2{n),  (11) 

where  the  computation  of  MaO  takes  time  proportional  to  T2(n). 

In  the  following  sections  we  show  how  the  quantities  N (),  Mi  ()  and  M2()  can  be  computed 
efficiently,  for  various  choices  of  Z. 


4  Conditioning  on  the  stock  price  at  selected  points. 


For  the  purpose  of  illustration  we  consider  first  an  approximation  based  on  conditioning  only  on  Hn, 
the  number  of  heads  (i.e.,  up-ticks)  up  to  time  n.  In  other  words  we  consider  the  approximation  (5) 
without  the  random  variable  Z.  Accordingly  we  suppress  the  argument  z  in  N  () ,  M i  ()  and  M2  () . 
It  is  easy  to  see  that 


where,  as  in  the  rest  of  the  paper,  we  interpret  the  binomial  coefficient  to  be  0  if  t  <  0  or 
A  <  0  or  f  <  h.  Each  of  these  expressions  can  be  computed  in  time  proportional  to  n.  Thus  from 
expression  (8)  we  see  that  the  time  required  to  compute  this  approximation  is  proportional  to  n*. 
Also,  from  (11),  the  error  bound  can  be  computed  in  time  proportional  to  n®. 

We  now  present  a  more  general  approximation,  where  we  choose  Z  to  be  the  combination  of  d 
random  variables 

for  some  chosen  values  of  wi  <  n2  <  ...  <  n.i  <  n.  Thus  the  approximation  for  an  Asian  call  is 
the  right  hand  side  of: 


E(A„  -  L)+  >  E 


(E(.4„| 


(12) 
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It  should  be  clear  that  we  get  increasingly  good  approximations  as  we  use  larger  values  of  d.  Of 
course  this  comes  at  the  expense  of  an  increased  computation  time,  as  we  see  below.  The  size 
of  the  sets  C{k),  in  expression  (6),  is  proportional  to  f(n)  =  n^.  As  mentioned  before  we  only 
need  to  show  how  to  compute  the  quantities  N  ()  and  M\  ()  efficiently.  Since  the  Z  in  this  case  is 
multi-dimensional,  we  write  N (fc,  hi, . . . ,  kj)  to  denote  the  number  of  paths  with  ff,,  =  k  and 

The  notation  Mi  {k,  hi ,  h2, . . .  ,  h  t,  h)  is  defined  similarly.  It  is  straightforward  that 


(Recall  that  we  are  interpreting  the  value  of  a  binomial  coefficient  (*)  to  be  0  if  h  <  0,  fc  <  0  or 
k  <  h.)  Also, 

Ml  (h,  hi ,  h2 , . . .  ,  hci;  t,  h)  = 

(N (h,  hi,  h2, . . . ,  he,  h,  he+i, . . . ,  hj)  if  n-e  <  t  <  We+i  for  some  1  <  e  <  d 
N{k,h,hu...,hu)  (14) 

'V(h,hi,h2,...  ,hi, h)  ifnd<t, 

which  represents  the  number  of  paths  through  (ni,  hi),  (712,  h2), . . . ,  [n.^,  hj),  (n,  k)  and  (t,  h). 
The  computation  of  this  expression  takes  time  proportional  to  Ti  (n)  =  dn.  From  expression  (8), 
the  total  time  required  to  compute  this  approximation  is  proportional  to 

T?f(n)Ti[n)  =  dn^^\ 

The  quantity  M2 (h,  hi ,  h2, . . .  ,  h f  1 ,  hi ;  (3,  h2)  can  be  computed  by  an  expression  similar 
to  the  one  above,  in  time  proportional  to  dn.  Thus  from  (11),  the  error  bound  can  be  computed  in 
time  proportional  to  dn4+®.  We  omit  the  tedious  details. 

5  Conditioning  on  the  average  of  the  underlying  random  walk. 

We  now  consider  the  approximation  (1)  where  we  take  Z  to  be  the  random  variable  Z,,,  which  is 
sum  of  the  random- walk  positions  from  time  0  to  time  n.  Thus  for  an  Asian  call,  our  approximation 
is  given  by  the  right  hand  side  of 

E(.4„  -  L)+  >  E  E(.4„  I  ff„,  Z„)  -L^ . 

Note  that  the  value  of  Zn  determines  the  average  position  of  the  underlying  random  walk  (not  the 
stock  price)  between  time  0  and  time  n.  This  in  turn  determines  the  geometric  average  of  the  stock 
prices,  which  is  simply  C,,  =  As  mentioned  in  the  introduction,  an  intuitive  reason  why 

conditioning  on  {H,,,  Zn)  might  yield  a  good  approximation  is  this:  paths  with  the  same  geometric 
stock-price  average  that  end  at  the  same  stock  price  will  likely  have  arithmetic  stock-price  averages 
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n 

Hn 

z„ 

Paths 

Gn 

Min.  An 

E(4,,|Z,.,7f„) 

Max.  An 

var(4„|Z„,  Hn) 

8 

6 

18 

4 

115.191 

115.575 

115.832 

116.346 

0.09895 

4 

0 

8 

100.000 

100.111 

100.167 

100.334 

0.00929 

10 

8 

33 

5 

120.893 

121.554 

121.907 

122.614 

0.16226 

5 

1 

20 

100.577 

100.666 

100.725 

100.812 

0.00404 

12 

10 

52 

125.978 

126.983 

127.416 

128.283 

0.22334 

6 

0 

100.000 

100.077 

100.150 

100.488 

0.00683 

Figure  1:  Illustrating  the  small  variance  of  the  arithmetic  stock-price  average  .4„  over  the  set  of 
paths  that  have  the  same  (a)  geometric  stock-price  average  (7„  and  (b)  up-ticks  .  The  parameters 
assumed  are  Sa  =  100,  L  =  90,  a  =  .2,  and  /j,  =  .05.  “Paths”  indicates  the  number  of  paths 
with  the  indicated  value  of  Hn  and  Z„.  “Min.  A„"  “Max.  .4,,,”  and  E(.4„|Z„,  if,,)  respectively 
indicate  the  smallest,  largest  and  average  value  of  .4,,  over  this  set  of  paths;  all  these  paths  have  the 
same  value  of  G,,.  var(.4„|Z„,  if,,)  is  the  variance  of  .4„  over  this  set  of  paths.  For  each  n,  the  first 
line  corresponds  to  the  values  of  Zn  and  fi„  that  give  rise  to  the  maximum  variance  of  .4,, ;  in  the 
second  line  we  choose  ii„  =  n/2,  and  so  that  the  number  of  paths  in  the  set  is  maximized. 


that  do  not  differ  too  much.  Figure  1  illustrates  this  point.  In  [10],  Rogers  and  Shi  consider  an 
approximation  based  on  conditioning  on  Zn  alone  rather  than  both  Z,,  and  ff„.  As  we  show  below, 
conditioning  on  both  random  variables  has  the  advantage  that  it  allows  efficient  computation  in  the 
binomial  tree  model.  Moreover,  the  resulting  approximation  is  strictly  better  than  the  one  obtained 
by  conditioning  on  Z„  alone. 

To  describe  the  computation  of  the  quantities  'V(),  Mi()  and  M2{)  we  introduce  some  useful 
notation.  Let  Ci{k)  denote  the  set  of  possible  values  of  Z,  on  paths  such  that  Hi  =  h.  In  other 
words,  Ci{k)  is  the  set  of  possible  values  of  the  random-walk  sum  Z,  on  paths  that  go  through  the 
lattice  point  (i,  k).  Note  that  G„(A;)  =  C{k).  Also,  let  Ni(k,  z)  denote  the  number  of  paths  with 
Hi  =  k  and  Z,  =  z.  Note  that  'Y„(fc,  z)  =  N {k.  z) .  Since  the  maximum  possible  absolute  value  of 
Z;  is 

l-t-2-|-3-|-...  +  i  =  i(i'f  l)/2, 

and  Z,  is  always  an  integer,  there  cannot  be  more  than  i(i  -f  1)  possible  values  of  Z,  .  Therefore  the 
sum  of  the  sizes  of  the  sets  Ci{h)  for  /i  =  0, 1, ....  t  cannot  exceed  i{i  -t- 1).  In  particular,  the  size 
of  Cn{k)  =  C{k)  is  at  most  a  constant  times  /(»)  =  n^. 

We  now  show  that  the  quantities  N (),  Mi  ()  and  M2O  can  be  computed  efficiently. 

First  we  consider  the  quantities  z).  We  assume  that  we  maintain  a  three-dimensional 
nx  n  X  n(n  -|-  1)  array  to  store  the  values  z):  the  first  two  dimensions  correspond  to  i  and 
h  and  the  third  dimension  corresponds  to  the  z.  Note  that  for  i  >  0  all  paths  with  Hi  =  h  have 
E;  =  2h  —  i,  and  Hi-i  is  either  h  -  lot  h.  In  either  case, 

Z,_i  =  Zi  —  Yi  =  Zi  —  2h  i. 

Thus  Ni{h,  z)  satisfies  the  recurrence 

Ni{h.  z)  =  Ni-i{h  —\,z  —  2h +  !)-{■  Ni-i[h.  z  —  2/i  -b  1). 

Starting  with  -Yo{0, 0)  =  0,  aU  the  AT, (ft,  z)  can  be  computed  inductively  (or  equivalently,  by  dy¬ 
namic  programming)  using  this  recurrence.  Since  there  are  n{n  +  l)/2  points  (i,  ft)  in  the  binomial 
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lattice,  and  each  such  point  has  number  of  possible  Zi  values  proportional  to  the  computation  of 
all  the  A’i  (h,  z)  takes  time  porportional  to  n*.  In  the  implementation,  we  first  compute  aU  (ft  ,  z) 
values,  which  are  used  in  the  computation  of  the  Mi  quantities,  as  we  now  describe. 

We  now  turn  to  the  computation  of  the  Mi  (ft,  z;  t,  ft),  the  number  of  paths  through  (t,  ft),  (n ,  ft) 
that  have  =  z.  We  can  partition  these  paths  according  to  their  Zt  value.  The  number  of  such 
paths  with  Zt  =  sis  equal  to  Ni{h,  s)  times  X  where  X  is  the  number  of  paths  between  (t,  ft)  and 
(n,  ft)  such  that  Yi+i  +l'i+2  +  . .  .+y„  =  z-s.  Note  that  all  paths  through  (t.  ft)  haveT)  =  2ft -t, 
so  X  is  the  same  as  the  number  of  paths  between  (0, 0)  and  [n  —  t,k-  ft)  such  that 

Zn-t  =  z-s  —  [n  —  t)Yt  =  z  —  s  —  [n  -  t)(2h  —  t). 

Thus 

Mi{k,z;  t.h)  —  ^  Nt{h,s)Nn-t(k  —  h.z  —  s  —  {n  —  t)(2h  —  t)), 

is€C-'((A) 

which  can  be  computed  in  time  proportional  to  Ti  (n)  =  v?,  since  the  size  of  C'i(ft)  is  no  more  than 
some  constant  times  Thus  the  total  time  required  to  compute  this  approximation  is  porportional 
to  n*  (for  pre-computing  the  Ni{h,  z)  values)  plus  expression  (8),  which  is 

n^f{n)Ti{n)  =  vJ . 

Since  n’  dominates  n*,  the  time  required  is  proportional  to  rtJ . 

We  mention  also  that  the  quantity  M^ik,  z:  (i ,  fti ;  t2 ,  ftj)  can  be  computed  using  similar  ideas, 
intimer2(n)  =  Thus  the  error  bound  can  be  computed  in  time  proportional  ton^/(»)r2(n)  = 


6  A  hybrid  algorithm 

We  now  briefly  describe  an  approximation  based  on  conditioning  on  a  random  variable  Z  that  is  a 
combination  of  the  two  random  variables  described  in  the  previous  two  sections.  That  is,  we  take 
Z  to  be  the  {d  -I-  l)-dimensional  random  variable 

Hni ,  ....  Hnj, 

This  gives  an  approximation  that  is  strictly  better  than  either  of  the  previous  two  approximations. 
While  the  details  of  the  computation  of  the  quantities  iV(),Mi()  and  M2()  are  rather  tedious,  we 
illustrate  the  basic  idea  by  means  of  the  case  d  =  1 .  That  is,  we  condition  on  the  random  variables 

Hn.Ht,Z„, 

where  Hs  a  fixed  integer  lying  strictly  between  0  and  n.  Our  approximation  for  an  Asian  call  is  the 
right  hand  side  of 

E(.4„  -  L)+  >  e[e(.4„  I  ftf„,  Ht,  Zn)  -L^. 

We  write  N (ft,  ft,  z)  to  denote  the  number  of  paths  with  =  z  that  go  through  (n,  ft)  and  (t,  ft). 
Also  we  write  Ml  (ft,  ft,  z;  ti,fti)  todenotethenumberof  paths  withJ?,,  =  z  that  go  through  (t,  ft). 
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(n,  k)  and  (ti ,  hi).  Recall  that  N (fc,  h.  z)  is  exactly  the  quantity  Mi  (fc,  t.  h)  that  we  computed 

in  Section  5  where  the  random  variable  Z  was  simply 

Now  we  describe  the  computation  of  Ml  We  assume  that  t  <  fi.  The  other 

cases  can  be  handled  similarly.  By  arguments  similar  to  those  used  to  compute  Mi  ()  in  the  previous 
section,  we  see  that 


Mi{k,h,z\  ti.hi)  = 

E  E 

s^Ct{h)  (/il) 


Nt{h,s)  X  (hi  —  h,  Si  —  s  —  (ti  —  t)(2h  —  t))  x 


Nn-ti  {k  -hi,  Z-Si-  {n-  ti){2hi  -  h)) 


,  (15) 


which  can  be  computed  in  time  proportional  to  n*  since  each  of  the  sets  Ct{h)  and  Cij(hi)  has 
size  at  most  some  constant  times  n^.  For  general  d,  the  computation  of  Mi  ()  takes  time  at  most  a 
constant  times  Ti  (re)  =  re^('^+*l.  Also,  the  size  of  the  sets  C(k)  in  expression  (6)  is  at  most  a  con¬ 
stant  times  /(re)  =  71*^+^ .  So  from  expression  (8),  the  time  required  to  compute  the  approximation 
is  proportional  to 

n^f(n)Ti(n)  = 

The  computation  of  the  quantity  M2()  can  be  carried  out  using  similar  ideas. 


7  Numerical  experiments 

Throughout  we  assume  that  the  option  expires  in  time  T  =  1,  the  number  of  steps  re  =  10,  and 
the  initial  stock  price  Sa  =  100.  We  vary  the  strike  L,  the  volatility  a,  and  the  drift  The  results 
are  displayed  in  Figures  2  and  3.  In  both  tables,  each  parameter  value  is  displayed  only  when  it  is 
different  from  the  one  in  the  previous  line.  For  ease  of  comparison  we  use  parameter  values  similar 
to  those  in  Rogers  and  Shi  [10]. 

Figure  2  shows  a  comparison  of  some  of  our  approximations  with  the  exact  value  of  an  Asian 
call  in  the  binomial  tree  model.  Fo  denotes  the  approximation  of  section  4  with  d  =  0,  i.e.,  with  con¬ 
ditioning  on  just  H,,.  Fa  denotes  the  approximation  of  that  section  with  d  =  3.  More  specifically, 
we  condition  on 

where  we  choose  re,  =  i[n/4j,  i  =  1, 2, 3.  Fo„  denotes  the  approximation  of  Section  5,  i.e.,  we 
condition  on  the  pair  [Zn,  Hn).  Finally,  F  denotes  the  exact  value  of  the  option  in  the  binomial 
model  given  by  E(.4„  -  L) + .  In  all  cases  we  see  that  the  agreement  between  F  and  Fa„  is  excellent. 

In  Figure  3  we  show  some  comparisons  of  our  approximations  with  the  Monte  Carlo  results  of 
Levy  and  Thmbull  [9].  Again  we  see  that  the  agreement  with  Vav  is  very  good. 
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Strike  L 

Vol  cr 

Drift 

n 

P 

ApxVb 

Apx  F3 

Apx  Vau 

Exact  V 

95 

0.05 

0.05 

1.0159 

0.6546 

7.1769 

7.1781 

i.mi 

7.1783 

100 

2.5941 

2.6860 

2.6939 

2.6939 

105 

0.1637 

0.2440 

0.2515 

0.2515 

95 

0.7819 

8.8151 

8.8151 

8.8151 

4.2617 

4.2859 

4.2881 

4.2881 

105 

0.7365 

0.7452 

95 

0.9739 

11.1115 

11.1115 

11.1115 

11.1115 

6.8080 

6.8080 

6.8080 

6.8080 

105 

2.5091 

2.5179 

2.5184 

2.5184 

1.0321 

0.5713 

11.9375 

11.9445 

11.9451 

11.9451 

3.3319 

3.5937 

3.6159 

3.6159 

110 

0.1485 

0.2462 

0.2552 

0.2552 

90 

0.09 

0.6350 

13.3859 

13.3883 

13.3885 

13.3885 

100 

4.6674 

4.8616 

4.8783 

4.8783 

no 

0.3227 

0.4726 

0.4853 

0.4853 

90 

0.15 

0.7310 

15.4151 

15.4154 

15.4154 

15.4154 

100 

6.8971 

6.9825 

6.9901 

6.9901 

110 

0.8727 

1.0934 

1.1098 

1.1098 

90 

0.50 

0.05 

1.1713 

0.4763 

16.1753 

17.2525 

17.3243 

17.3252 

100 

10.8140 

12.1138 

12.2208 

12.2232 

110 

7.2794 

8.3115 

8.4067 

8.4082 

0.4890 

17.0228 

18.0257 

18.0943 

11.5780 

12.8308 

12.9343 

12.9365 

7.9038 

8.9217 

9.0173 

90 

0.15 

0.5081 

18.2776 

19.1687 

19.2301 

19.2308 

100 

12.7439 

13.9163 

14.0138 

14.0158 

110 

8.8794 

9.8660 

9.9577 

9.9591 

Figure  2;  Comparison  of  approximations  with  the  exact  value  of  an  Asian  call  in  the  binomial  tree 
model 


Strike  L 

Vol  a 

Drifts 

V3 

Monte  Carlo 

95 

0.10 

0.09 

8.86 

8.90 

8.91 

8.91 

100 

4.66 

4.86 

4.89 

4.91 

105 

1.75 

1.94 

1.96 

2.06 

90 

0.30 

14.45 

14.90 

14.93 

14.96 

100 

7.93 

8.70 

8.77 

8.81 

no 

3.77 

4.54 

4.60 

4.68 

90 

0.50 

17.02 

18.03 

18.09 

18.14 

100 

11.58 

12.83 

12.93 

12.98 

no 

7.90 

8.92 

9.02 

9.10 

Figure  3:  Comparison  with  Monte  Carlo  results. 
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